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Abstract 



The survival probability of large rapidity gaps in pp collisions is calculated 
for several different eikonal models of the Gaussian form. Results obtained 
for models based on partonic interactions are quite similar. The Regge-pole 
model predicts a higher value of < l^ 2 > . 



1 Introduction 



It has recently been suggested |], that the observation of a large rapidity gap in the 
r] — 4> lego plot, constructed for exceedingly high energy p-p interactions, may serve as a 
signature for W-W fusion associated with the production of a Higgs boson. The practical 
utilization of this idea as a useful trigger for rare electroweak processes depends on ones 
ability to reliably assess the survival probability < \S\ 2 > . This is defined as the 
fraction of events for which spectator events do not fill the rapidity gap of interest. 

The survival probability is easily defined in the eikonal model in impact parameter 
space. We use amplitudes normalised so that 

f = -|/(M)| 2 (1) 
a tot =47r/m/(s,0) (2) 
a(s,b) = ±- fdqe-^ h f(s,t) (3) 
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From which we derive the b-space formulae [[| : 

Otot = 2 J dblma(s, b) (4) 

a el = J db\a(s,b)\ 2 (5) 

s-channel unitarity implies that \a(s, b)\ < 1, and when written in a diagonalized form we 
have 

2Ima(s,b) = \a(s,b)\ 2 + G in (s,b) (6) 
from which we obtain for the inelastic cross section 

a in = J dbG in (s,b) (7) 

s-channel unitarity is most easily enforced in the eikonal approach where, assuming that 
a(s,b) is purely imaginary, we can write 

a(s,b) =i(l - e - Qis ' b) ) (8) 

where the eikonal Q(s, b) is a real function. 

Our assumption that a(s,b) is purely imaginary is not compatible with analyticity and 
crossing symmetry. These are, easily restored upon substituting s a s a e~ l ~ . From 
Eqn. (7) we can express Gi n (s, b) as a function of Q(s, b) 

G in (s,b) = l-e-™^ (9) 
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Note that P(s,b) = e 2n ( s ' b ) is the probability that no inelastic interaction takes place at 
impact parameter b. 

We follow Bjorken 0] and define < \S\ 2 > as the normalized multiplication of two 
quantities. The first is a convolution over the parton densities of the two interacting pro- 
jectiles presenting the cross section for the hard parton-parton collision under discussion. 
The second P(s,b) is the probabilty that no other interaction takes place in the rapidity 
interval of interest. In the eikonal formalism we have: 



, / qh(s, b)P(s, b)d 2 b 



<\S\ 2 >= J rfMh (10) 
J clh{s, b)d A b 

where clh(s, b) denotes the amplitude associated with hard collisions that can be expressed 
through the eikonal Q,n{s,b) using Eqn. (8). Some preliminary calculations of < \S\ 2 > 
have been presented in Ref . |2j . Following this pioneering effort there were also a number 
of attempts made to estimate < \S\ 2 > using Monte Carlo techniques ||, ||. It is important 
to note that these assessments of the survival probability are model dependent. A number 
of models are available [6 - 11] which provide a good reproduction of the data in the ISR 
- Tevatron range. As we shall show these models differ in their estimates of fi#(s, b) and 
P(s,b) in the high energy limit of LHC and SSC. It is therefore pertinent to carefully check 
the dependence of < \S\ 2 > on the phenomenological input required in Eqn. (10). This is 
the main aim of this note, where we have attempted a systematic study of < |5| 2 > and 
it's sensitivity to the input parameters. 

In the following we assume that both an{s,b) and Q(s,b) in P(s,b) are well approxi- 
mated by Gaussians 

a H (s,b) = v H (s)e ^ (11) 

2Q(s,b) = 2u(s)e~^ (12) 

This input assumption has been verified by the analysis of Ref. fl(| [111 , where it has been 
shown that eikonal models of this form provide an excellent reproduction of the cross 
section data and particle distributions in the energy range 5 < y/s < 1800 GeV. The 
main advantage of assuming the input in a Gaussian form is that the integration in Eqns. 
(4) and (10) can be carried out analytically, whence the total cross section 

-l) n -V(s) 



db(l-e- n ^) = 2nR 2 (s)Y / 



nln 



l 

= 2nR 2 (s)[lniy(s) +C- Ei(-v(s))] (13) 

Ei(x) denotes the integral exponential function Ei(x) = J*^ ~dt , and C is the Euler 
constant ( C = 0.5773 ). For v(s) > 1 we have 

a tot -> 2nR 2 (s) [lnv{s) +C- e^ (s) ] (14) 
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The inelastic cross section is given by 

a in = J db(l - e - 2n(s - b) ) = uR 2 {s)[ln(2u(s)) +C- Ei{-2v(s))\ (15) 
and therefore a e \ = a to t — Qin = 



= nR 2 {s)[ln{^) + C + Ei(-2u(s)) - 2Ei{-v{s))) (16) 
Again for u(s) > 1 we have 

a el -> TiR\s)[ln(^-) + C + e - 2u{s) - 2e~ u{s) ] 

One can extract the values for R 2 (s) and u(s) using the expressions for a tot and a d 
given in Eqns. (13) and (15). As expected at the Tevatron energy where the cross 
sections are known, the spread in values of R 2 {s) and v(s) for the different models is 
rather small. It is only when the model parameters are extrapolated to higher energies 
does the difference become significant, and allows one to test the theory on which the 
parametrization is based. We shall estimate the importantance of this dependence when 
making predictions for LHC and SSC energies. 

The integration of Eqn. (10) yields: 

<\S\ 2 >=a{^Y 1 {aM (17) 
2v 

where a — and 7(0, x) denotes the incomplete gamma function 

H 

~f(a,x) = f ' z a - x e~ z dz (18) 
Jo 

For u(s) ^> 1 Eqn. (17) simplifies to [0 : 

, aT(a) a „ t s 

u{s) a u{s) 

A general numerical mapping of < |*5| 2 > as a function of a and v is shown in Figs. 1 
and 2. 

A realistic asssessment of the survival probability is subject to a considerable ambi- 
guity, due mainly to the following reasons: 

1) There is no clear definition of the hard component of the b-space scattering amplitude 
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or eikonal. Moreover, it is not clear whether the growth of atot with s is due to the soft 
or to the hard sector. Bjorken has suggested that one can estimate R 2 H from the low 
energy cross section where a tot (pp) rs 40 mb. Hence R 2 H is energy independent and a > 
1. Implicit in this estimate is the assumption that a to t growth is associated with the soft 
sector. This point of view is not common to all models. In particular, models based on 
parton interactions [§, [UJ associate the growth of a to t with the hard or semi-hard sector. 
Hence R 2 H is energy dependent and a ~ 1. Bjorken |J based his numerical estimates on 
the minijet model 0. Indeed in this type of model R 2 H approaches a constant in the 
high energy limit. Here, the growth in a to t is due to gluon-gluon interactions which are 
semi-hard. Having no better insight to this problem, we follow the recipe suggested by 
Bjorken, and fix R 2 H from the low energy data at ^/s ~ 10 GeV. 

2) The input information required in Eqn. (17), i.e. a and 2v is obtained from fits to the 
rising pp and pp cross sections. This rise is a consequence of both a process of blackening, 
i.e. a rise of v(s) , and expansion, i.e. an increase of R 2 . However, these two processes 
compensate each other as is readily seen in Eqn. (14). The dependence of < \S\ 2 > on 
these parameters is completely different as is evident from Eqn. (19). As we have noted, 
there are a number of phenomenological models which reproduce the available data well 
in the 10 < -/s < 1800 GeV energy range. These models differ in their estimates of a 
and 2z/, and we wish to examine the stability of < |*S| 2 > with respect to variations of the 
input parameters a and 2v , as deduced from the phenomological models. 

3) We would like to emphasis that all our estimates of < \S\ 2 > are based on the eikonal 
approach, i.e. on the assumption that the representation for the scattering amplitude 
given in Eqn. (8), is a valid approximation at high energies. This form allows us to 
formulate the expression for the survival probability in a simple and transparent manner. 

In the following we elaborate on these points by discussing a few input models: Our 
results are summarized in Table I, where we list the values for i/, R 2 H , R 2 , atot and < IS') 2 > 
for energy values of the Tevatron, LHC and SSC. 

1.1 Minijet-parton model 

In this model || the eikonal for the partonic i-j collision is given by 



n( S ,b)(s,b) = W ij (b)a„ 



.QCD 



(20) 



where a- 



~ s 



j-i 



. For Wij , the Chou-Yang formulation is assumed W2 



/4 



(fiijb) 3 ''K 3 (nijb) 



(21) 



96tt 
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Eqn. (21) is reasonably well approximated by a Gaussian. To realize this we have fitted 
Eqn. (21) numerically, utilizing the parameters for fcj given in Ref. || . To estimate 
the hard component, we follow the method suggested by Bjorken and fix R 2 H from the 
model predictions at yfs = 10 GeV, where Otot ~ 40 mb. This is a natural choice as in 
this model, we have relatively small changes of R 2 (s) with energy. Not suprisingly, our 
estimates of < |5"| 2 > are compatible with those of Refs. 0. f|. 



1.2 Regge-pole model 

An impressive reproduction of the experimental total cross section data is obtained @ 
by utilizing a simple Regge-pole model. For our high energy analysis we are interested in 
the super-critical Pomeron where the amplitude is given by: 

f(s, t) = iCe B % t s a ®- 1 sin[0.5ira(t)) (22) 

with a(t) = 1 + e + a't . Donnachie and Landshoff fit the data with values of the 
parameters C = 21.7 mb and e = 0.0808. In addition we utilize a global fit || to B, the 
nuclear slope, and obtain R% =5.2 GeV~ 2 and a' = 0.2 GeV" 2 . The b-space transform 
of Eqn. (20) is 

a(s, b) = ^C-^exp[-^-] • [J2?«n(£a(0) + Z) - ?-a>cos£a(0) + Z)\ (23) 



2\(3\ 2 4/3 



where 



R\ = Rl + a'lns 

7r 4 r/ 2 

P\ 2 =Rt + 
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na'b 2 , , 

Z = — i ^ 24 

8\(3\ 2 v ; 

As both e and a'Cl, Eqn. (21) is well approximated by a Gaussian ( see Eqn. (11) ) 
with 

The model in the form suggested by Donnachie and Landshoff is not appropriate for 
fitting data at high energies, as for very high energy it violates unitarity (a(s, b = 0) > 1) 



5 



above y/s ~ 5 TeV. For our evaluation we use a very similar eikonalized version suggested 
by Cudell and Margolis @, with C = 24 mb, e = 0.093 and a' = 0.25 GeV~ 2 . For the 
hard sector our choice is less straightforward than before, as R 2 (s) ~ In s, and there is no 
obvious way of defining the hard component. In Table I we present R 2 H as the threshold 
of Eqn. (25). Note that the results obtained from the Regge-pole model are in marked 
contrast to the other models investigated in this note. The difference will be elaborated 
apon later. 



1.3 Lipatov-like Pomeron 



A simple parametrization for the Lipatov-like Pomeron |L3[], has been suggested in Ref. 
rP~CTl - The eikonal is given by 



"<•■*> = (s£- e " TO < 26 > 

where the following two parametrizations of R 2 (s) provide good fits to the data: 

^ = 04 + 05(^5)°" (27) 

R 2 L2 — a 4 + a 5 Vlns + a§lns (28) 



l L2 

dj are fitted parameters. In this model a regular Pomeron with trajectory ct(0) = 1, 
is appended to the Lipatov-like Pomeron. Again the choice of what to take for R 2 H is 
ambigous, as R 2 (s) is energy dependent and it's low energy limit is exceedingly small. We 
adopt an arbitrary definition as suggested in Ref. P] and use R 2 H to be the value at y/s 
= 10 GeV, where a to t ~ 40 mb. 



1.4 Dual parton model 



This is a multi-component model describing soft and semihard multiparticle produc- 
tion. The eikonal is given by 

Q(s, b) = Q s (s, b) + Q H (s, b) - Q TP (s, b) - Q L {s, b) (29) 

where the last two terms correspond to the triple Pomeron and loop contributions. As 
f2rp(s,b) and f^(s,b) are reasonably small , this is effectively a two component model 
whose parameters are: 
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where 



OS 



Cs e 



R 2 S = B + a' Ins 

with C = 40.8 mb, e = 0.076 and a' = 0.24 GeV~ 2 



(31) 



oh 
8nR 2 H 



e iR n 



(32) 



where R 2 H 



B. In this type of model R 2 H is defined as the low energy threshold limit of 
R 2 s (s) given by Eqns. (30) and (31). The an is calculated in lowest order QCD and is 
dependent on the value taken for the p™ m cutoff. Numerical values for an at different 



energies are given in [11 . 



2 Conclusions 



Our results are summarized in Table I. The results obtained for the various partonic 
models fllP"! , JTTJ] for < \S\ 2 > in the LHC-SSC energy range are remarkably stable. 



We note that the models of |T0[ and |TTJ , even though very different in their construction, 
yield rather similar input parameters, as summarized in Table I. Ref. || differs from 
the above, in that it has the highest values for v. These are compensated for by having 
the corresponding lowest values for a, producing final results which are similar to those 
obtained in Refs. 



10 and 



The survival probability obtained from the Regge-pole model [8-10] are considerably 
higher. On examining the input parameters used in our calculation, we find that the 
difference can be traced to the fact that the v values associated with the Regge-pole 
model are the smallest. Hence, in order to be compatible with the data the model requires 
relatively large values of R 2 H and R 2 which give rise to high < \S\ 2 > . 



We conclude with a more general comment. Clearly the questionable aspect of a 
calculation such as presented here, is the fact that the definition of the hard component 
is not unique. We have followed Bjorken's suggestion [fj], and have fixed R 2 H from the low 
energy data. It is likely that our estimates for R 2 H are on the conservative side, so that in 
reality one could expect even higher values for the survival probability than are given in 
Table I. 
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Table I Parameters and predictions of different models 



TV K 11 

Model 




i/(s) 


D2 

K H 


r>1 

R A 


Cfot 


< \s\ 2 > 




TeV 




Gey 2 


GeV z 


mb 


% 


Mmijet (6) 


1.8 


2.50 


14.41 


19.74 


72 


13.2 




16.0 


3.90 


14.41 


22.60 


107 


5.5 




40.0 


4.75 


14.41 


23.20 


121 


3.8 


(1-7 Q \ 

Regge [< yj 


1.8 


1.11 


25.41 


35.80 


76 


32.6 




16.0 


1.48 


25.41 


40.16 


102 


22.1 




40.0 


1.68 


25.41 


41.99 


117 


18.1 


Lipatov f ^ 


1.8 


1.60 


15.78 


25.39 


75 


19.7 




1U.U 


9 fiQ 
z.uy 


lO. ( o 


98 n9 


1 1 1 
no 


R 9 




40.0 


3.44 


15.78 


29.37 


134 


4.9 


Lipatov 2 ^ 1U ) 


1.8 


1.44 


16.19 


28.64 


76 


20.6 




16.0 


2.24 


16.19 


32.66 


115 


9.2 




40.0 


2.77 


16.19 


34.23 


137 


5.8 


Dual parton ( n ) 


1.8 


1.83 


10.56 


28.47 


75 


9.6 


KMRS[B-2] 


16.0 


2.23 


10.56 


32.67 


109 


5.3 




40.0 


2.43 


10.56 


34.43 


124 


4.0 
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Figure captions 

Figure 1: Contours of percentage of the survival probability < IS*] 2 > as a function of 

a and v . 

Figure 2: Graph of log ( % survival probability ) versus a for selected values of v . 
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